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Abstract 

The capture of electrons by the nucleus Be from the three-body initial state p + e~ + 7 Be in the 
continuum is studied. On the basis of the expansion of the three-body continuum wave function in 

1 /2 

a small parameter e ~ (m e /m p ) ' [m e (m p ) is the electron (proton) mass], the role of the protons 
on the electron capture is considered. The results are compared with the traditional treatment of 
the electron capture by the nucleus 7 Be. For stars with the density and temperature like in the 
center of the Sun the studied mechanism can make non-negligible contribution to the capture rate. 
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I. INTRODUCTION 



The process of the electron capture by the nucleus 7 Be is important since it contributes to 
the low energy part of the spectrum of neutrinos radiated by the Sun. Besides, it is obvious 
that the balance of the disappearance channels of 7 Be in the Sun regulates the amount of the 
nucleus 8 B which is the source of the high energy solar neutrinos. This is the main reason 
why this process attracted considerable attention over many years Q, [3, [3, 3, 6, 7]. 
Practically all the discussion so far of the electron capture in 7 Be is concentrated on consid- 
ering the electron wave function in the vicinity of the nucleus and on the screening effects 

Q 

on it; the study of this capture in the plasma was done in Ref. [7J]. 

In what follows we will estimate the role of the process that is usually not included in the 
standard theory of the pp cycle in the Sun. Let us first note that in the standard theory of 
this cycle the destruction of the nucleus 7 Be takes place in the following binary reactions 

p + 7 Be -> 8 B + 7, (1.1) 
e~ + 7 Be -> 7 Li + v. (1.2) 

Since the nucleus 7 Be participates in both processes, instead of the binary reactions (11. ip 
and (11. 2p we consider the contribution to the electron capture rate from the three-particle 
initial state p + e~+ 7 Be. In this case the following reactions can take place 

/ 7 Li + v + p, (1.3) 
p + e~ + 7 Be -> 8 B + 7 + e" , (1.4) 
\ 8 B + e~. (1.5) 

As it was shown in Ref. |2J the screening corrections for the electrons in the continuum are 
rather small. Therefore, we consider in the initial state the bare Coulomb interaction in all 
two-body subsystems e~ +p, p+ 7 Be, and e~+ 7 Be. In this case one can immediately realize 
that there is a qualitative difference between the binary and ternary mechanisms* of the 
electron capture. Indeed, if one starts from the three-body initial state, then the processes 
(II. ip and (11.21) should be interdependent because the wave function of three charged particles 
cannot be presented as a product of pair wave functions, as is required by the binary processes 



* We use the terms binary- and ternary reactions as synonyms for the reactions in the two- and three 
component systems. 
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(11. lp and (jl.2p ' . Here we apply the method and results obtained in Ref. 9] where an 
alternative to the Born-Oppenheimer approach has been suggested for the calculations of 
wave functions in the continuum for three charged particles. 



II. CALCULATIONS 

In analogy with the approach developed in Ref. [9J , the continuum wave function of three 
charged particles can be expanded in a small parameter e, 

l(M + m p )(m p + m e ) 

where, in addition to the electron and proton mass, also M, the mass of the nucleus 7 Be, 
enters. The expansion of the wave function of the studied three-body system is then 

V(r,R) « ^ {r,R) + e^x{f,R) + ••• . (2.2) 

Since the parameter e w 0.0233 one expects the effects of the second term at the right hand 
side of Eq. (12. 2 p at the level of 2 % in comparison with the first term. 
It was found in [s| that in the limit e — > the Jacobi coordinates r and R (see Fig.[T]) 
separate. It means that the structure of the wave function ty (r,R) is 

tt„(r, R) = V C (R) * c (f, Z = Z x + Z 2 ) , (2.3) 

where ty c (R) is the Coulomb wave function describing the relative motion of the proton 
and the nucleus 7 Be, and *f? c (r, Z = Z\ + Z 2 ) is the Coulomb wave function that describes 
the motion of the electron in the field of the effective Coulomb potential of the charge 
Z = Zi + Z 2 . The crucial point is that the wave function \J/ C7 (r, Z = Z\ + Z 2 ) depends on 
the distance between the electron and the center of mass of the subsystem of heavy particles. 
It means that even if the distance between the electron and the nucleus 7 Be is zero, as is 
required by the Hamiltonian of the weak interaction, the wave function ^ (r, Z = Z± + Z 2 ), 
defining the probability of the electron capture by the nucleus 7 Be, should be taken at a 
non-zero distance |r| = /3\R\, where (3 = 1/7 is the ratio of the proton and 7 Be masses. 

T There is only one exception corresponding to the case when all three particles are at very large distances 
between themselves Q, which is not applicable to the electron capture. 



(2.1) 
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FIG. 1: The Jacobi coordinates for the p + e + Be system. 



It is clear that this phenomenon appears due to the electron movement in the Coulomb field 
of two charged particles with positive charges. Following the arguments presented above 
we now consider two effects acting in the opposite directions. On one hand, increasing the 
effective positive charge of the heavy particles system by one unit will enlarge the electron 
capture rate. On the other hand, using the Coulomb wave function at finite distances instead 
of the function taken at zero distance should damp the capture rate. 

Taking into account that the nuclear matrix elements of the reactions (11.211 and (11. 3B are the 
same, as a measure of influence of the third particle (the proton in this case) on the capture 
rate of electrons by the nuclei 7 Be we introduce the ratio q(R,T), which is the function of 
the distance R between the particles and the temperature T, 

c(RT] = f oo mPR,Z = 5)\*e- E / kT dE 

q{ ' ' / °° F(Z = 4, v) e~ E / kT dE ' 1 ] 

Here the denominator on the right-hand side of the equation contains the quantity that 
enters the electron capture rate from the continuum for the reaction (jl.2p [1]. The Fermi 
function F(Z, v) is given by the equation 

F(Z,u) = 2nu/(e 27TU -I), (2.5) 

where the parameter v is given by the equation v = —Zam e /p = —Za/v, a is the fine 



structure constant, and p = ^2m e E is the electron momentum. It is obtained by using the 



io|. 



solution of the Dirac equation with the Coulomb potential 
Analogous integral in the numerator should reflect the effect of the Coulomb potential on 
the electron in the continuum for the reaction ( 11. 3ft . For the wave function ^(r, Z), we use 
the Coulomb continuum wave function for the state with zero angular momentum 

*S(P) = (2-6) 



where the function F (r],p) satisfies the equation 

d 2 F 2 V 



+ [1- -L]F = 0, (2.7) 



dp 2 p 



with p = p/3R, and 77 = —Zaym e /E = —Za\f2jv is the Sommerfeld parameter. The func- 
tion Fo(r],p) can be expressed in terms of the Kummer function M (see Ref. |ll|], Chap. 14) 

as 

F (v, p) = C ( V ) p e-v M(l - i v , 2, 2%p) , (2.8) 

where C^(r]) = 2nr]/(e 2 ^ - 1). 

Instead of the quantity q(R,T), one can consider 

' J / °° |^g(0,Z = 4)| 2 e- s / fcT rfE ' 1 J 

Using Eq. (12.71) . one obtains 

\^ C E (0,Z)\ 2 = 27r V /(e 2 ^-l), (2.10) 
which is of the same form as the Fermi function (12. 5ft ; however r\ = \[2v. 



One can find the following integral representation (see Ref. Chap. 13) 

M(l-irj,2,2ip) = ^ /' e 2 ^f— V* dt , (2.11) 

TIT] Jo \ t J 

for this function. In the numerical calculations of the integral over energy in Eqs. (12.41) and 
(12. 9p . we used this representation of the Kummer function for the energies E > 0.1 keV. 
In the interval E < 0.1 keV, the function under the integral in Eq. (l2.1ip strongly oscillates 
which makes the calculations difficult. Instead, we applied the program PFQ developed in 
Ref. [3]. Let us note that for the energies E > 0.1 keV the program PFQ and Eq. (12.1 lj) 
provide the same results to a high degree of accuracy. 
We also introduce the mean value < ^(-Ro, T) > of the function s(R, T) 

<q(Ro,T) >= N J e M o q (R, T) dR = 4irR*N J e a x 2 s(RqX, T) dx , (2.12) 



where 

j 2 

5? _ ^"-"-o 



(ZzMl . 4vrjR 3 1 



iV- 1 = / e "^dR = -^J-[l + v / 27i(l - 0( (2.13) 

j ye y2 

is the normalization constant. In Eq. (12.131) . the function <f)(y) is the error function (see 
Ref. [llj, Chap. 7). It can be seen from Eq. (12.121) that the quantity < s(Rq,T) > depends 



on the mean distance Rq between the particles defined by the density in the Sun and on the 
temperature T. The mean value < sc{Ra,T) > is defined analogously by using qc{R,T). 
We checked the precision of numerical calculations of the quantities < s(R ,T) > and 
< ^c(Ro,T) > by using Mathematical and also independent numerical procedures. The 
results of these two independent ways of calculations agree within the required accuracy, 
which is 0.1 %. 

III. RESULTS AND DISCUSSION 

The results of the calculations are presented in Table I and Figs. [2]- [H In Fig.[2l the depen- 
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FIG. 2: The dependence of the mean value < s(Rq,T) > on the temperature T and the value of 
the mean distance Rq. Solid line, Ro = 0.1 x 10 4 fm; dashed line, Rq = 0.25 x 10 4 fm; dotted line, 
Ro = 0.5 x 10 4 fm; dashed and dotted line, Rq = 1.0 x 10 4 fm; dashed and double dotted line, 
Ro = 1.5 x 10 4 fm. 

dence of the mean value < s(Ro, T) > given in Eq. (I2.12p on the temperature T and the value 
of the mean distance R is shown. A weak dependence of < q(R ,T) > on the temperature 
means, as it follows from Eq. (12 .4p . that the temperature dependence of the electron capture 
by 7 Be is almost the same for the ternary and binary reactions. Such a behaviour can be 
understood from the fact that in both cases only rarely is all the kinetic energy carried by 
one particle, which is the electron. On the other hand, the dependence of < q(Ro,T) > on 
the value of Rq shows that the contribution to the capture rate of the ternary reaction is 
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presumably suppressed in stars but it can be at the same level as the contribution to the 
capture rate for the binary reaction or even prevail over it at very high densities. This is 
natural because at short distances between the particles the factor of the larger effective 
charge acting on the electron will dominate. The same conclusion can be drawn from Fig.0 
In this figure, the solid and dashed curves practically coincide. This again shows a very 
smooth dependence of < <;(Rq, T) > on the temperature. Let us note that the values of 
the i?o, considered in Fig. [2], correspond to rather dense stars. For example, the value of 
Ro = 10 4 fm corresponds to the proton density p p = 1673 g/cm 3 , which is about 11 times 
larger than in the center of the Sun. Let us further discuss the electron capture by 7 Be 
solely in the Sun. 




-4 
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FIG. 3: The dependence of the mean value < s(Rq,T) > on the values of the mean distance Rq. 
Solid curve, kT = 1.161 keV, (3 = 1/7, and Z = 5; dashed and dotted curve, kT = 1.161 keV, 
j3 = 4/7, and Z = 6; dashed curve, kT = 1.5 keV, (3 = 1/7, and Z = 5. The dashed and dotted 
curve corresponds to analogous calculations for the nuclei 4 He. 

In Table I, we show the influence of protons on the electron capture in the Sun in more 
detail. For the Standard Solar Model, we choose Model SSMBP2004 [13]. According to 
Fig. 6.1 of Ref. jljj], the maximal intensity of the electron capture by the nuclei 7 Be takes 
place at the distance R S /R Q ~ 0.06, where R Q is the radius of the Sun, and it drops to 
one half at R s /R & ~ 0.03 and R s /R & ~ 0.1. Using the data on the temperature, the 
density, and the fraction of the hydrogen in this area of the Sun, we obtain the mean values 
< (;(Rq,T) > and < $c(Ro,T) > presented in Table I. In the second column, we add the 
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mean value calculated at R s /R & = 0.007 which is close to the very centre of the Sun. It can 
be seen that the change of the average quantities < s(Ro,T) > and < sc{Ro,T) > is very 
smooth. 

Table I. The mean values < <;(Ro,T) > and < <;c{Ro,T) > for the electron capture by the 
nuclei 7 Be in the Sun. 



<s(Ro,T) > 


0.0991 


0.0991 


0.0965 


0.0913 


<q c {Ro,T) > 


0.0718 


0.0717 


0.0696 


0.0658 


R s /R & 


0.007 


0.03 


0.06 


0.1 


kT (keV) 


1.353 


1.300 


1.161 


1.088 


P P (g/cm 3 ) 


52.1 


51.6 


48.9 


44.9 


R x 10~ 4 (fm) 


3.179 


3.188 


3.250 


3.340 



As can be seen from the first row of Table I and from Fig.Hl the contribution to the capture 
rate of the ternary reaction at the Sun is about 10 % of the binary one*. This means that 
it should increase sensibly the burning out of the nuclei 7 Be in comparison with the binary 
reaction, thus decreasing the concentration of the nuclei 8 B that appear after the capture of 
protons by 7 Be. 

Comparing the first and the second rows of Table I shows a difference of 3 % between the 
values of < s(Ro, T) > and < <Jc(-Ro 5 T) >. This variation arises from the difference between 
the relativistic and nonrelativistic estimations of the electron wave function at zero distance 
for the binary reaction (11. 21) . 

In Fig. H] we show the variation of < ^(Rq, T) > for the reaction ( 11. 3ft for larger intervals of 
T and R . 

One can consider in analogy the influence of the nuclei 4 He on the electron capture by 7 Be 
in the Sun. In this case, Z — 6, /3 — 4/7, and at the radius R s /R & = 0.06 the mean 
distance between the nuclei 4 He is Rq = 5.34 x 10 4 fm. Then one obtains from Eq. (I2.12p 
that < <;(Rq,T) >=0.0036, which is about 27 times smaller that the analogous value of 
< ^(i?cb T) > for the protons given in the fourth column of Table I. Evidently, this influence 
on the electron capture is negligible. This can also be seen from Fig.O 
The main conclusion following from our calculations is that the three-body process due to 
the presence of the proton in the vicinity of the nucleus 7 Be results in the capture of the 

* The next term in the expansion of Eq. (|2.2p is expected to change it only by w 2 %. 
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FIG. 4: The dependence of the mean value < $(Rq,T) > on the temperature T and the value of 
the mean distance Rq. Solid line, i?o = 3.25 x 10 4 fm; dashed line, Rq = 2.75 x 10 4 fm; dashed and 
dotted line, R = 3.75 x 10 4 fm. 

electron by an effective charge Z = 5 instead of Z — 4, which is qualitatively the new effect 
that cannot be simulated by introducing the Debye screening. This effect can increase the 
rate of the electron capture from the continuum by 7 Be in the Sun, which will reduce the 
concentration of the nuclei 8 B that appear after the capture of protons by 7 Be. 
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